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1. Introduction 

The large iV problem of usual gauge theory remains a formidable problem, despite the 
existence of a closed set of large N equations, the well-known Makeenko-Migdal equations 
0. Some progress has been made recently by Migdal @, although it seems that addi- 
tional input is needed in order to finally solve these equations. Interestingly enough, these 
equations can be solved in two dimensions, as shown long ago by Kazakov and Rostov 
||. Unfortunately, as soon as one introduces dynamical scalar particles or quarks in the 
adjoint representation, the large iV problem again becomes intractable. For some refer- 
ences on this topic, see || . In light of recent progress in four dimensional supersymmetric 
Yang-Mills theories, initiated in the work of Seiberg-Witten || and furthered in || and 
0, it is tempting to ask whether supersymmetry helps in solving the large iV problem. 
Our goal in this paper is a modest one, instead of working in four dimensions, we ask 
whether supersymmetry helps in solving the large N problem in two dimensions, where 
introducing dynamical adjoint matter already complicates the problem a lot. The answer 
is yes, though the method adopted here is completely different from that of |J. The fact 
that our model is the dimensional reduction of a three dimensional N = 1 super-gauge 
theory or even of a four dimensional super-gauge theory may hint at possible simplification 
of the large N problem in these models. 

If one starts with deriving an equation in the super-gauge model parallel to the or- 
dinary Makeenko-Migdal equation, one need to derive more equations in order to get a 
closed system. One soon realizes that infinitely many equations are needed, so this way of 
proceeding is hopeless. As it turns out, the only equations we need to supplement the MM 
equations are the ones resulting from the Ward identities associated with supersymmetry. 
These identities are valid only when supersymmetry is not dynamically broken. This will 
be demonstrated in sect. 4, where we put the model in a spatial box. 

In the pure gauge theory, as being solved in [0 , the only relevant modes are topological. 
The solution of Wilson loops with intersections is quite nontrivial. In addition to the usual 
area law, the dependence of these loops on areas of windows is polynomial without a definite 
sign. This implies that if one tries to formulate any string theory (as attempted at in ||), 
one would have to introduce fermions on the world sheet. Indeed a formulation of such 
theory has proven quite unwieldy. It may appear surprising that the solution of Wilson 
loops in the supersymmetric theory is simpler than that in the pure gauge theory, as will be 
seen in sect. 3. It appears that the planar Wilson loops are described by the Nambu string 
without folds. It remains to see whether at the string loop level (1/N corrections) the 
correspondence persists. In any case, our result already indicates the following interesting 
picture. In two dimensions, when there is only gauge field, namely the theory is purely 
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bosonic, then one has to introduce fermionic degrees on the world sheet. However in the 
super-gauge theory, where there is a fermion in spacetime, one has only bosonic degrees of 
freedom on the world sheet (the fold- less constraint can be easily implemented). Further 
study is necessary to understand other loop-like quantities in addition to the Wilson loop. 

To begin with, let us write down the action of the N = 1 supersymmetric Yang-Mills 
theory in which the super- multiplet consists of a gauge field A M , an adjoint scalar 0, and 
an adjoint real fermion A, each field is a Hermitian matrix. Here for simplicity we consider 
a U(N) gauge group which makes no difference than a gauge group SU(N) in the large N 
limit. We shall follow conventions in ||. Let a = a = — 1, a 1 = —a 1 one of the Pauli 
matrices. The supersymmetric action is 

S = j d 2 xtr [~Fl v ~ \{D^) 2 - iXa^D^X - gXa 3 [^ X]j , (1.1) 

where g is the coupling constant. The field strength is defined according to F^ u = d^A v — 
d v A^ + ig[Afj,,A v ] and the covariant derivative of a Hermitian matrix field A (in ( |1 . 1|) it 
is either <fi or A) is defined by D^A = d^A + ig[A^, A}. The action ( |lTT|) is invariant under 
the following supersymmetry transformation 

SAp = -2zAa M e, 
5<p = 2zA(T 3 e, (1.2) 
5X = a x F Q1 e - a^a 3 D^e. 

For large iV considerations, it is often convenient to rescale all fields such that the 
action is weighted by a factor N, also we need to hold g 2 N fixed for large N. Thus let 
gVN -> g, and A^ -> (l/g)A fM -> (y/N/g)A IM , X -> (VN/g)X and (/> -> (VN/g)(f). Since 
all fields are rescaled by the same factor, the transformation law in ( |1.2|) remains the same, 
while the action is now weighted by a overall factor N/g 2 : 

S=^J d 2 xtv ^-±F 2 u -±(D^) 2 -iX^D^X-Xa 3 [<P,X}y (1.3) 

In the definition of the field strength and the covariant derivatives there is no explicit 
dependence on g. Now it is g 2 not g 2 N held fixed in the limit iV — > oo. 
We will work in Minkowski spacetime throughout this paper. 

The plan for the rest of the paper is as follows. In sect. 2 we shall consider a set 
of Makeenko-Migdal equations, and a Ward identity associated to supersymmetry. This 
Ward identity, together with one of Makeenko-Migdal equations, does not yet form a closed 
set of equations for the Wilson loop and a quantity with two insertions of the fermion 
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field and a third quantity. The validity of the Ward identity depends upon unbroken 
supersymmetry, which we will prove in sect. 4. We then proceed in sect. 3 to argue that the 
third quantity is indeed vanishing, so we have a closed set of equations. These equations 
are easily solved. Sect. 4 can be ignored if the reader does not wish to read the proof of 
unbroken SUSY. A special gauge is chosen in sect. 4 to discuss the Hamiltonian formulation 
of the theory, in order to keep a supersymmetry on a spatial lattice. For gauge group 
SU(N), the Witten index is calculated to be tr (-1) F = 1 or tr (-1) F = (-1)^. The 
ambiguity in determining the sign of the index is discussed and resolved. In either case, it is 
nonvanishing and signaling unbroken SUSY. Sect. 5 is devoted to a discussion. The model 
studied in this paper is shown to be a dimensional reduction of N = 1 three dimensional 
super-gauge theory in appendix A, where we also show that N = 2 super-gauge theory 
in two dimensions is a dimensional reduction of the N = 1 four dimensional super-gauge 
theory. Another set of MM equation and Ward identity is discussed in appendix B. 



2. Equations of Motion of Large N Wilson Loops 

As usual the Wilson line associated to a curve C xy with end points at x and y is 
defined by 



U{C xy ) = Pexp ^ jT A^dx^j , 



which transforms under A^ -> UA^U' 1 + id^UU' 1 as U(C xy ) -> U(x)U(C xy )U~ 1 (y). 
For a closed loop C xx , W(C XX ) = -^tr U{C XX ) is gauge invariant. Its expectation value is 
what we want to calculate. 

Another gauge invariant quantity relevant to our discussion is obtained by inserting 
the fermion field A at two points x, y on the loop C. These two points divide the loop into 
two segments of curves C xy and C' yx . Now W\(C xy , C' yx ) as a matrix is defined according 
to 

{W x (C xy ,C' yx )) a p = itr X a (x)U(C xy )X^y)U(C yx ), 

W\ is a two by two matrix. 

One's first instinct is to write down the usual Makeenko-Migdal equation derived from 
the identity 
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[dAJ0dA]tr U(C xx )e lS = 0. 



The equation is presented in appendix B. It is easy to see that this equation, unlike the 
MM equation in the pure gauge theory, will involve three different quantities. To get a 
closed set of equations, more Schwinger-Dyson equations are needed. Proceeding further, 
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one will soon realize that the number of equations will never terminate, namely a closed 
set of Schwinger-Dyson equations will involve infinitely many equations. 

Supersymmetry plays an important role in this model. Clearly, if SUSY is not broken 
dynamically, there are many Ward identities one can write down. We shall prove in sect. 4 
that indeed SUSY is not broken, therefore for an arbitrary functional of fields F(A, <p, A) 



J [dAd<j>d\}5 e {F{A, </>, \))e iS = 0, (2.1) 



where S e (F(A, <fi, A)) is the SUSY transformation of F. The reason behind the above iden- 
tity is the following. First of all, the action is invariant under supersymmetry transforma- 
tion, then the measure is invariant tooS. If the vacuum is annihilated by the super-charges, 
then the boundary conditions in the path integral is also invariant under SUSY, thus 

[cL4 e # e dA e ]F(A e , \ e )e iS = J [dAd(f)dX]F(A, (f>, X)e lS 

implying ( |2.1|) . If SUSY were broken, it would still be possible to write certain Ward iden- 
tities. One would include effects of non-invariance of boundary conditions, such identities 
obtained may be called "anomalous" Ward identities. 

Let F = tr X(x)U (C xx ) in ( |2.1[ ), a Ward identity is readily written down, using trans- 
formation law ( |1 . 2j ) : 

(tr F^UiC^a 1 -^ D^{x)U(C xx ))a»a 3 + 2 I dy^(tv \(x)U{C xy )\{y)U{C' yx ))a^ = 0, 

(2.2) 

where the last quantity is what we have already introduced. This is an equation of two by 
two matrix. Next, use the relation []TT[ 



and the relation 



tr U(C XX ) = iti F^U{C XX ) 



itr D^UiC^) = d^tr <j>{x)U{C xx ) = d^W^C^) 



the Ward identity ( |2.2| ) is written as 

" ^^(C-)^a 3 + 2 j dy^W x (C xy , C'^a^ = 0, (2.3) 



1 Without explicit calculation, the measure is invariant at least up the the first order in e, since 
the Jacobian is bosonic and SUSY transformation is linear. 
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where we use d/da(x) to denote d / da 01 (x) , there is only one independent area element in 
two dimensions. The above is an equation relating three different quantities including the 
Wilson loop. To solve this equation, we expand the matrix W\ as follows 

W x (C xy , C' yx ) = + W 2 ia 2 + W 3 a 3 . (2.4) 

Under Lorentz rotation A — > exp(^cr 1 )A, cr^ transforms as a vector, a 2 and a 3 transform 
as a scalar. To keep Lorentz invariance, one then demands W 7 ^ transform as a vector, W 2 
and Ws transform as a scalar. Under parity reflection A — > a 3 A, a 3 is a scalar and a 2 is 



a pseudo-scalar. Thus, W3 is a scalar and W 2 is a pseudo-scalar. Substituting (|2.4|) into 
(p73|), one deduces 

q^) W ( C ™) = -2ifdy^W»(C xy ,C' yx ), (2.5) 
dpW^C^) = 2 j> (e^dy»W 2 + ^ v dy v W z ) , (2.6) 
= 0, (2.7) 



where the anti-symmetric tensor e^ u is specified by eoi = 1. These equations are valid 
even for a finite N. The first two equations tell us that in order to calculate W and 
it is enough to know W\. The last equation says that the total flux of TU M along the loop 
is zero. This is important for us, it allows us to extend the quantity 

$(C xy , C' yx ) = I" dtpW^Cg, C' yx ) (2.8) 

as a function of x and y into inside the loop C, by deforming the contour as shown below. 





Figure 1 

Note that this extension of $ as a function of y depends on the location of x. The 
above argument is bit hand- waving. A more formal argument is the following. For a simple 
loop, a loop without intersection, $ as a function of y when x fixed, is well-defined on the 
loop C. It is then always possible to analytically extend it onto the interior C. Again, 
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the extension depends on the position x. For a loop with intersection points, $ may be 
multi-valued when y is at one intersection. The reason is that when it crosses around the 
loop, it will cross this intersection point at least twice. Formula (|2.7|) does not guarantee 
$ be unique there. However, as we shall see, that $ satisfies a differential equation which 
is not multi-valued anywhere when x itself is not an intersection point, and its solution 
can not be multi-valued anywhere. We thus believe that <3> is single valued even at an 
intersection point y, as long x is not that point. 

From definition ( ^.8| ), it follows that W M = d y ^<&. $ is well-defined globally, it follows 
that OqWi — d\Wo = 0, consistent with eq.(|2~7|). Stokes formula when applied to ( |2~5p 
yields 

g^ W ( C ™) = 2i j da(y)d y »d y ^, (2.9) 

where the area integral extends to the domain enclosed by the loop C. Now it is desirable 
to derive an equation for $. 

Instead of deriving the ordinary MM equation associated to translational invar iance 
in the gauge field, we start with an equation associated to translational invariance in the 
fermionic field A, 

J [dAd<pd\]tv j^U(C xy )\(y)U(C yx y s = 0. (2.10) 

The derivative when acts on X(y) gives rise to a delta function 5 a p5 2 (x — y) and a product 
tr U(C xy )tr U(C' yx ). Each factor in the product is not gauge invariant unless x coincides 
with y. This is guaranteed by the delta function factor. Two additional terms result from 
the action of the derivative on exp(zS'): 

2N 

d 2 (x - y)(tr U(C xy )tT U(C' yx )) + — <r"0 M (tr \(x)U(C xy )X(y)U(C' yx )) 

2lN a 3 (tr [<t>,\](x)\(y)U(C>)) = 0. 



9 2 

This equation of two by two matrix is valid for an arbitrary N. In the large iV limit, apply 
the factorization theorem to the first term 

5 2 (x - y)W(C xy )W(C yx ) + ^a^Wx - ^a 3 W^ x = 0, (2.11) 



where the new quantity 



W^ x = l(tr [<P,\}(x)U(C xy )\(y)U(C yx )). 
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To make use of the matrix equation ( |2.11| ), we make the expansion, as in (j2.4|) 

W^x = + W 2 i<J 2 + W 3 a 3 . (2.12) 



Substituting this expansion and that in ( |2.4| ) into the matrix equation ( |2.11| ) and reading 
off coefficients of each basis matrix, we obtain 



d » Wfx = _ z ^ 3 + g ls 2 {x - y)W{C xy )W{C' yx ) 
e^d tl W„=tW 2 , 
dxWa - d W 2 = iWi, 
diW 2 -d W 3 = -iW . 



(2.13) 
(2.14) 
(2.15) 
(2.16) 



As we showed earlier, the vector W M is curl-less, so we conclude from (|2.14| ) that W 2 = 0. 
This is a pseudo-scalar. Eqs.( pTT5|) and (|2.16| ) do not respect Lorentz invariance unless 
d^,W 2 = 0. This implies that W 2 = const. W 2 is also a pseudo-scalar, so if it is independent 
of positions x and y, the only reasonable constant is zero. Thus, W 2 = W 2 = 0. As far 
as Wfj, is concerned, there is still an unknown quantity W3 in ( |2.13| ) . If one can show that 
this quantity is also vanishing, then eq. Q2.13D together with the Ward identity ( |2.5[ ) forms 
a closed system of equations for and W. 



3. Solution of Wilson Loops 

To make eq. (p~5|) or ( |2.9| ) together with ( |2.13| ) a closed set of equations, the central 
problem is to determine W3, a scalar quantity. It was already pointed out in the previous 
section that the pseudo-scalar W 2 = 0. 



It is seen from the expansion (|2.12|) that a non- vanishing W 2 would have measured the 
disparity between the two off-diagonal elements of W^x- Its vanishing says that there is no 
disparity. Similarly, a non-vanishing W3 measures the disparity between the two diagonal 
elements of W^x- It is natural to guess W3 = 0. An exchange between the diagonal 
elements can be achieved by transformation 

A — > ct 1 A, (3.1) 



which exchanges Ai and A2. Indeed this is a discrete symmetry of our theory (|1.3| ), provided 
that a simultaneous transformation </> — > —<p is made. This is because a 3 changes sign under 

(PD- 
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Now we can draw strong results from this symmetry. First consider W\, it should 
be invariant under the discrete symmetry (3.1). However, both a 2 and a 3 change sign 



under this transformation. It follows from expansion (2.4) that W2 = W3 = 0. We 



already argued for W2 = in the previous section with the help of eqs.( |2.15D and ( |2.16| ). 



Substitute W2 = W3 = into those equations, we find = 0. This result can be obtained 
by observing the quantity W^x too. Under the discrete transformation, W^x changes its 
sign, however do not change their sign, therefore they must be zero. Unfortunately, we 
can not conclude W2 = W3 = from this symmetry, for a 2 and cr 3 do change their sign. 
Nevertheless, as a consequence of SUSY Ward identity and eq.( |2.14| ), W2 = 0. Thus, we 
have inferred that all coefficients in expansion of W^x are zero except W3. 

We make the conjecture that W3 = 0, which is very natural in our opinion. The fact 
that W2 = does not follow from any symmetry of the model encourages us to make this 
conjecture. It is plausible that W3 = is related to W2 = by certain duality. The latter 
is a pseudo-scalar while the former is a scalar, duality usually relates a scalar quantity to 
a pseudo-scalar quantity. For example, the electro-magnetic duality relates the vector Ei 
to the pseudo- vector Bi. We shall discuss another set of MM equation and Ward identity 
in appendix B, where we present an argument which is close to a proof of W3 = 0. 

It is possible that both W2 and W3 become singular when points x and y all approach 
an intersection point. But such complication will not alter our result obtained below, as 
long as we stay away from intersection points. 

With W3 = and = eq.( [2.13|) together with eq.( [2.9|) forms a simple system 



of equations 

-^W{C xx ) = 2i I da(y)d y »d y ^, 



8v»K* = - v 52 ( x - y)W(C xy )W(C' yx ) 



(3.2) 



2 

Now it is a simple matter to solve the Wilson loop from the above equations. One simply 
substitutes the second equation into the first one, and performs the area integral. If the 
loop C is smooth at x and not an intersection point, half of contribution of the delta 
function is picked up, because the area integral is restricted inside the loop, one then has 

' -W{C xx ) = -4-W{C xx ). (3.3) 



da(x) v 2 

The loop C may have many intersection points, therefore many windows. It is reason- 
able to assume that W{C) depends on the loop only through areas of these windows. This 
is indeed dictated by (|3.3|). If the point x is on a segment of the loop separating a window 
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Si from the infinite area outside the loop, a variation of Sa(x) is simply a variation of area 
Si. Eq. (|3.3|) results in 

d Si W{S) = --^W{S). (3.4) 

If the point x is sitting on a segment separating two windows Si and Sj, and Sj is inside 
Si, then 

(d j -d i )W(S) = --^W(S). (3.5) 

From the second equation in (|3.2|) , $(x) is solved as $ ~ ln(a; — y) 2 W(C xx ), as 
long as x not an intersection point. $ is singular at y = x. This singularity can be 
regularized as usual by introducing the ie prescription. By definition, = d y v<& ~ 
(y — x)fj,/(y — x) 2 W(C xx ). This result can be contrasted with a perturbative consideration. 
Without coupling to the gauge field and the scalar, the first factor agrees with the usual 
Dirac propagator. The factorization when it is coupled to bosonic fields is interesting. We 
like to caution ourselves that this result may not be valid at an intersection point x = y. 

vanishes at a non-intersection point y = x, if the ie prescription is used. This fact will 
be used in the discussion in appendix B. 

Back to Wilson loops. Eqs.( |3.4[ ) and ( |3.5| ) are enough to determine the function W(S). 
Let us consider a few examples. The simplest one is a simple loop without intersection 

■ 2 

point, as shown below. The solution to (^J) i s W(S) ~ exp(— ^-S). The proportional 
coefficient must be one, for when S shrinks to a point W = 1. Note that the area law 
is exactly the same as in the pure gauge theory [|J, although here we work in Minkowski 
spacetime. 




Figure 2 

The next example is the 8 shaped curve. The two windows are not separated by a 
segment, we need use only ( |3.4|) . The result is 

W(Si, S 2 ) = exp (~{Si + S 2 )) • 

Again this is an area law agreeing with the pure gauge theory. So far both Wilson loops 
respect the Nambu string behavior. 
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Figure 3 

A crucial example is the loop in the next figure, where £2 borders on Si. Eq. (|3.4|) 
gives diW = —i(g 2 /2)W, this together with (|3?5| ) gives d 2 W = —ig 2 W. Thus the solution 
to these equations is 

W(S 1: S 2 ) = exp ^—(S! + 2S 2 )) = exp (-^(Si + S 2 + S 2 f) . 

There is no power dependence on S 2 , unlike in the pure gauge theory. The above formula 
is perfectly in accordance with the Nambu string. 




Figure 4 

The final example is shown in the figure below. There are three windows. Our 
equations then determine the Wilson loop 



in*, . ,s 2 . ,%) = exp I -\(Si + 2S 2 + 3S S ] 



exp ( — 2~(^i + S 2 + S 3 + S 2 + S 3 + S3] 



also agrees with the Nambu string. 
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Figure 5 

It is not hard to convince oneself that the standard area law should persists for all 
kinds of loops. It is therefore very interesting to learn that in a theory of more complicated 
spacetime physics, the world sheet picture is simpler. Certainly we are not claiming here 
that the whole theory is described by the Nambu string without folds, since there are many 
other physical operators independent of Wilson loops. 



4. Hamiltonian Formalism and the Witten Index 

We discuss the Hamiltonian formalism in this section for two purposes. First, we 
want to determine whether supersymmetry is dynamically broken. To adapt an argument 



of Witten in [|I0[ to calculate to Witten index, we need to put the system into a spatial 
box of finite length. Second, to study nonperturbative effects, it is often tempting to put a 
system on a lattice. If one wants to make use of supersymmetry, this can not be normally 
done on a full spacetime lattice, since supersymmetry implies translation invariance in 
two directions. However, some models can be put on a spatial lattice without spoiling 
a subset of supersymmetry generators, provided no spatial translation is generated by 
these generators. This can be done in two dimensions for N = 1 supersymmetry. In four 
dimensions, one requires at least N = 2 |L2[ . 



For a gauge theory, there is one more complication. One need to fix a gauge in the 
Hamiltonian formalism. A gauge must be chosen such that it is invariant under some of 
supersymmetry transformations. In the two dimensional N = 1 Yang-Mills theory at hand, 
there are two supersymmetry generators Q a , each is Hermitian and satisfies Q 2 a = 2if, 
H is the Hamiltonian. One can not keep both Q a , since the anti-commutator of the two 
generators gives rise to the spatial translation generator. Now if one chooses the temporal 
gauge Aq = 0, it is easy to see from the transformation law ( |1.2| ) that no supersymmetry 
survives this gauge. However 

5(A + </>)= 2i\(a s - l)e = -4zA 2 e 2 , 
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so if 62 = 0, the combination Aq + <f> is invariant under Q\. We thus fix a gauge in which 
Aq + (j) = 0. With this gauge choice 

S=^ Jd 2 xtr (^(doAf + doAD^+^do^ + iXdoX + iX^D^ + Xil-a 3 )^^}^ 

9 (4-1) 
where A is the spatial component of the gauge field, and D x is the spatial covariant 
derivative. The Yukawa coupling involves only A2 in the above action. 

The canonical momenta are read off from the action 

^ (d A + D x 4>) , 

9 N < 4 - 2 ) 
IIa = i^X. 

The Hamiltonian is then 

H=^ Jdxti Q(9 A) 2 + ^(d <p) 2 - iXa 1 D x X + A(a 3 - !)[</>, A]^ . (4.3) 
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It is also straightforward to write down the unbroken super-charge Q 



1 



Q x = 2 J dxtr + A 2 (n A - ^D x 4>)^j . (4.4) 

To check the relation Q 2 = 2H, one should notice the fact that since n>, is the same as A, 
the anti-commutator is the half the value of the usual anti-commutator. Specifically, 

We will not try to write down the Hamiltonian in terms of the link variable and other 
fields, except making a comment on the role of fermions. There will be a doubling problem 
as usual on a spatial lattice. Here one solves the problem by putting Ai on even sites, and 
A2 on odd sites, much like what is done in fll"2"|| . In the super-charge (|4.4j), although 



will be only assigned on even sites, the term D x <p involves (p both at an even site and an 
odd site, it is easy to check that the relation Q\ = 2H is satisfied by this prescription of 
solving the doubling problem. The continuum limit is naturally achieved. 

In the remaining part of this section, we calculate the Witten index. We follow closely 
a calculation by Witten of the index in the four dimensional super-gauge theory in ||1U|| . 
We put the system into a spatial box with boundary x = 0, L, with periodic boundary 
conditions. If the Witten index is nonvanishing for all finite L, it is certainly nonvanishing 
in the infinite volume. One may consider a gauge group U(N). But since all fields are 
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in adjoint representation, the U(l) sector is free and then does not affect the issue of 
supersymmetry breaking. So we will consider gauge group SU(N) (If there is a matter 
sector, the U(l) sector can not be ignored, and indeed is a subtle problem as investigated 
in ||10|| .) Unlike in four dimensions, it seems that the weak coupling does not come to our 
rescue, the reason is that g 2 has a dimension of mass squared, therefore the meaning of 
weak coupling is senseless. Still, we have a dimensionless combination g 2 L 2 . When L held 
fixed, we can make this combination arbitrarily small. If the Witten index tr (— 1) F is not 
zero for a small g 2 L 2 , then it is not zero for arbitrary g 2 L 2 since the index is invariant 
under arbitrary deformation of parameters. 

Now the trick of Witten consists in reducing the problem to a problem of zero modes 
(with vanishing momentum). Any nonzero momentum mode will have a energy greater 
than 1/L. For a zero mode, we will see that an excited state associated to the zero modes 
of A has a energy g 2 L, which is much smaller than 1/L if g 2 L 2 is small enough. Thus it 
is safe to ignore nonzero momentum modes in the discussion. Write the Hamiltonian in 
terms of canonical momenta of fields 



H=jj dxtr I -ITA + -Ui) I ( J.r.tr U A D, 



(4.5) 



+ ^ J dxtr (^(D x 4>) 2 -iXD x X + X(a 3 -l)[^X} 

The last term of ( f4.5| ) tells us that the zero modes satisfy 

D x <j> = 0, D X X + i(a 3 - 1) [<f>, A] = 0. (4.6) 

We still have freedom to do spatial gauge transformation, and it is always possible to gauge 
transform A into a constant matrix -4.(0). This constant matrix can not be gauged away 
in general, since the Wilson line tr exp(zA(0)L) can not be gauged away with periodic 
gauge transformation. Because only the zero mode A is left, the second term in ( |4.5| ) is 
independent of non-zero momentum modes of <p. This is why we can focus our attention 
on the zero mode of <p i n the first place. Now the solution to the first equation in (fOJ) is 
given by 

(j)( x ) = e - iA ^ x (j)(0)e iA ^ x . 

The periodic boundary condition <p(L) = (f)(0) implies that [</>(0), exp(zA(0)L)] = 0. Thus, 
the Hermitian matrix and the unitary matrix exp(iA(0)L) can be simultaneously diago- 
nalized, which in turn implies [A(0), (fi(0)} = 0, therefore <p(x) = </>(0). With remaining 
constant gauge transformations, we can always put A(0) and <p(0) into a maximal Abelian 
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sub-algebra of SU(N). Let A(0) = E a A^ a , 0(0) = E a <M a - fa are generators of this 
sub- algebra. 

The second equation of (^4.6|) can be solved just as the first one, and the result is that 
A a (0) = J2 a Kyt a ■ Now apparently there is a vacuum which is annihilated by IT^ and 11^. 
Let it denoted by |fi). The fermion part form a Clifford algebra {A™, A^} = l/2S a pS a b. Its 
representation is 2 N ~ l dimensional. To construct this representation, let A a = A" — zA^, 
and A a = A" + Now {A a , X b } = {A a , A 6 } = and {A a , A 6 } = 5 ab . We let |fi) also be 
the vacuum annihilated by A a , then any other state can be written as A" 1 . . . A ai 

To count the true vacua, our final ingredient is the residual gauge group. First, one 
can do gauge transformation U = exp(i2irt a x/L). It is periodic and shifts A a by an amount 
of 2tv/L. So A a is a periodic variable. While 4> a is not restricted. Second, there is global 
gauge group consisting of global gauge transformations mapping the maximal Abelian sub- 
algebra into itself and is called the Weyl group. For SU(N), it is the permutation group 
Sjv- Any physical vacuum is invariant under a permutation. Since II A a = ~ *l/^A a an< ^ 
A a is periodic with a period 2n/L, an excited state, according to (fOD, has an energy 
g 2 L. As long as g 2 L is much smaller than 1/L, non-zero momentum modes can be safely 
ignored. 

The problem is that there are excited states of arbitrarily small energy, for </> a 's are 
not restricted. One may put a cut-off on the space of 4> a , say tr </> 2 (0) = ^ a 4> 2 < A 2 . This 
cut-off is gauge invariant and will cause no trouble to have the theory well-defined. Now 
let A to be large enough such that g 2 / A 2 is much smaller than 1 (so that the excited state 
associated to <p a has an energy g 2 /(A 2 L) much smaller than then the spectrum is 

discrete, the Witten index is well-defined. In the end of calculation, one can push A to 
infinity without changing the Witten index. 

One can either assume that \Q) is invariant under permutations, or pseudo-invariant 
(changes its sign under a odd permutation). If it is invariant, then no more invariant states 
can be constructed. This can be shown along the line of [ I0| . In this case tr (— 1) F = 1. 
If one assumes |0) be pseudo-invariant, an invariant state can be generated by acting on 
it by the pseudo-invariant operator A 1 . . . A^ -1 . This state , call it has the statistics 
(— 1)^, if one assumes that the pseudo-invariant state \Q) is fermionic. The Witten index 
is then tr (— 1) F = (—1)"^. This is an ambiguity hard to resolve in the four dimensional 
super-gauge theory jl0|1 . We argue that this issue can be resolved in our model. Note 
that if is invariant, then |f2) is pseudo-invariant and has a statistics (— l)^ -1 . This 
state is annihilated by all A, so as a state it should be treated on the equal footing as 

This means that when |0) is pseudo-invariant and |fi) is invariant, the latter should 
be regarded as a bosonic state, and the former has a statistics (— l)^ -1 not —1 as we 
previously assumed. So the Witten index is 1 in this case too. 
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We thus have shown that the Witten index is nonvanishing for a finite L and small 
g 2 L 2 , therefore it is nonvanishing for arbitrary L and g 2 . Supersymmetry is not broken in 
the 2D super-gauge theory. If one further demands the existence of vacuum in the limit 
N = oo and its statistics being well-defined, one has to choose tr (— 1) F = 1. The other 
choice ( — 1)^ does not make any sense. 

5. Discussion 

Tracing the reason why the equations of motion in the super-gauge theory are con- 
siderably simpler than those in the bosonic theory, we realize that the Schwinger-Dyson 
equation associated to translation of the fermion field is a first order differential equation, 
since the kinetic term of the fermion in the action is of the first order. The Ward identity 
( |2.3| ) relates the Wilson loop to the loop with two insertions of the fermion field in a simple 
manner. Recall that much difficulty in dealing with the ordinary MM equation arises from 
the second order derivative of the Wilson loop in the equation, which is absent in equations 
discussed in sect. 2. 

We have seen that the solution of the Wilson loop is really the Nambu string without 
folds. One would like to proceed further to study other physical observables, in order to 
learn more about the string theory underlying this super-gauge theory. Evidently, the 
string theory possesses a spacetime supersymmetry, and the physical spectrum should 
furnish a representation of SUSY. 

The power of combining the Schwinger-Dyson equations, in the guise of loop equations, 
and Ward identities associated to SUSY is manifest in our model. We have studied one set 
of these equations in the previous sections. In appendix B, we shall study another set of 
equations, where interesting results are also obtained. Although the method we present in 
this paper is markedly different from the holomorphy technique and duality argument of 
Seiberg- Witten, we suspect that there is intimate relationship. This may become evident 
if we study a high dimensional super-gauge theory, for only there duality also enters into 
loop variables The N = 1 2D super-gauge theory is a dimensional reduction of a 

N = 1 3D super-gauge theory (appendix A), thus with additional input, hopefully the 
large iV problem in this model can also be solved. We plan to study this model in the 
future. 

A more straightforward application of considerations here would be to the N = 2 2D 
super-gauge theory. As shown in appendix A, there is a complex scalar in this model, and 
one more Majorana spinor in the adjoint representation. Classically, there are many vacua, 
characterized by a moduli space, very similar to N = 2 four dimensional theories studied 
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in || . And the kind of "duality" suggested in sect. 3 becomes obvious in this model. Also, 
the N = 1 4D super-gauge theory dimensionally reduces to this model, one may learn 
things in the high dimensional model by studying this 2D model. 

Finally, we have studied only the Wilson loop and a couple of related physical observ- 
ables. Physical problems such as the spectrum in this model are still open. To solve them, 
one would need study more physical observables, a systematic scheme would be valuable. 
Indeed a possible such scheme, the free variable representation of master fields, has been 
the subject of a flurry of recent activities [H| - [FS|. The master field was constructed by 
Singer for the 2D pure gauge theory. The relative ease in constructing it is due to the free- 
ness of the master field in a special gauge. The model studied in this paper then presents a 
challenge: The master field is no longer expected to be free for different momentum modes. 
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Appendix A 

The N = 1 four dimensional super-gauge theory without matter contains a vector 
super- multiplet, in which there are a gauge field A m and its super-partner A. Here following 
conventions in || Latin letters m and n are used to denote the spacetime index. In the 
so-called Wess-Zumino gauge, there is an auxiliary field D, which is also in the adjoint 
representation. The action 

S = J d 4 xtv (-\Fl n - iXa m D m X + l -D 2 ^ (Al) 

is invariant under the SUSY transformation 

6 A m = -iXa m e + iea m \ 
S\ = a mn F mn e + iDe, (A2) 
5D = D m Xa m e + ea m D m X. 

A three dimensional action is readily obtained by dropping out x 2 as well as A2. 
To have a supersymmetric theory, one demands that A is a Hermitian matrix, instead a 
complex one. It is easy to see that this is consistent with transformation (A. 2), where e 
also becomes a real spinor and SA2 = 0. In addition, 5D = 0, since all a m except a 2 are 
symmetric. Thus, the auxiliary field can be dropped out. Now, our two dimensional super- 
gauge theory as given in ( |1 . 1| ) is simply a dimensional reduction of the three dimensional 
super-gauge theory. 

In going from four dimensions to three dimensions, we dropped out A2 and half of 
degrees of freedom in A. If one does not do so, a super-gauge theory can still be obtained, 
where A2 becomes a scalar in three dimensions, and A decomposes into two Majorana 
fermions, A = Ai + 1X2- This is a iV = 2 super-gauge theory in three dimensions, and 
the auxiliary field is also kept (5D 7^ 0). Reducing one more dimension x 3 , a A = 2 
super-gauge theory in two dimensions is obtained. The field content is: A gauge field, two 
Majorana spinors, one scalar <pi = A3 and one pseudo-scalar <p2 = A%. Since the action is 
a direct reduction of (A.l) into two dimensions, we will not write it down here. This is a 
model of great interest to study, along the line of this paper. The "duality" we mentioned 
in sect. 3 becomes apparent in this theory, the exchange of role of <pi and <p2- If one forms 
a complex scalar from these two field, then the duality transformation is just a complex 
conjugate transformation. Two spinor fields will also get exchanged too. Classically, there 
are many vacua corresponding to different expectation values of the complex scalar. So 
this model is similar to the N = 2 super-gauge theories in four dimensions. In addition, as 
we have seen, this is a dimensionally reduced model of N = 1 4D super-gauge theory. So 
if one wishes to probe some physics in this 4D theory, the N = 2 2D super-gauge theory 
should serve as a good starting point. 
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Appendix B 

A Schwinger-Dyson equation associated to translation of the fermion field is considered 
in the main body of this paper. It is not the original Makeenko-Migdal equation, which 
is associated to translation of the gauge field. It is easy to generalize the ordinary MM 
equation in the pure gauge theory to our model, which derives from 
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[dAd<i>d\}tr -^-^UiC^y 3 = 0. 



Taking terms of the fermion field and of the scalar field in the action into account 
(tr D v F^(x)U(C xx ))-i{tx [<P,D^]U(C xx ))-2(tv \a»\(x)U{C xx )) 

(B.2) 



+ g 2 N J dy»5 2 (x-y)(tr U(C xy )tr U(C' yx )) = 0, 



where the integral in the last term is taken as properly regularized when x = y: The 
contribution of the delta function is ignored except when x = y is an intersection point 
and both U(C xy ) and U(C yx ) are nontrivial. Applying the factorization theorem in the 
large N limit to the last term in (B.2) 

du da^(x) W{Cxx) + W ^ Cxx) + 2iW ^ C ™) + % 9 2 j dy^i* - y)W(C xy )W(C' yx ) = 0, 

(B.2) 

where 

W; = ±(tr [<j>,D^<j>]{x)U{C xx )) 

and W^(C XX ) is a term of W\ in the expansion (|2.4| ) and when C' yx = 0. In addition to 
the Wilson loop, the generalized MM equation involves two additional quantities. When 
x is not an intersection point, the last term does not contribute, for the integral is regu- 
larized @. Now W(C XX ) depends only on areas of windows of the loop C, the derivative 
d/da(x)W(C xx ) must be independent of x in the vicinity of x, as long as x is not an 
intersection point @. We deduce from the MM equation (B.2) the following identity 

W%(C XX ) + 2iW^C xx ) = (fl.3) 

This equation will play a crucial role in checking the consistency of our conjecture Ws = 
below. We expect that (B.3) is no longer true when x is an intersection point, otherwise 



2 If x is an intersection point, this derivative is discontinuous across x, since the area derivative 
on a different side involves different window areas. 
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the MM equation (B.2) would be identical to the one in the pure gauge theory, where the 
solution for W(C) is completely different from what we obtained in sect. 3. 

Next consider another Ward identity associated to supersymmetry. Substitute F = 
tr [(f), A] (x)U (C xx ) into the general identity (27T) and make use of the SUSY transformation 



( |1.2|) , we obtain a two by two matrix equation 

(tr [(P,F 01 ]{x)U(C xx ))a l - (tr [0, D^](i)[/(C M ))a' t a 3 - 4i(tr XX(x)U(C xx ))a 3 

(BA) 



+ 2 j (tr [(f>, \}(x)U{C xy )\{y)U(C' yx ))a^ = 0. 



All the terms except the first term are familiar quantities. The first term is vanishing, 
if x is not an intersection point. The reason is simple. Note that tr [(f>, Foi](x)U(C xx ) = 
tr (4>Fqi(x)U(C xx ) — (f>(x)U(C xx )F Q i(x)), and the effect of F i(x) when acting on U{C XX ) 
is an area derivative, whether it acts from the left or right, so this quantity is identically 
zero. After this term is dropped out from the above equation, the equation is rewritten in 
terms of quantities defined before 

2 £ dy»W^ x (C xy ,C yx )a^ = 4tW x (C xx )a 3 + W%(C xx )ay. (B.5) 

Remind ourselves that this equation is valid only when x is not an intersection point. Now 
we are in a position to check the conjecture W3 = 0. As argued in sect. 3, the nonvanishing 
terms in the expansion of W\ in terms of sigma matrices are W^. They are vanishing too, 
if the segment C' yx = and properly regularized, see sect. 3. From this fact and (B.3), it 
follows that W^(C XX ) = 0. We then see that the r.h.s. of (B.5) is zero, which in turn 
implies 

I dy» (W 2 ia 2 + W 3 a 3 ^ ^ = 0, 

where we used the expansion ( |2.12|) and the result = 0. This equation is consistent 
with W2 = W3 = 0. Indeed this is not even too much a weaker equation, as it might 
appear. It is equivalent to the following two equations 



(j)dy°W 2 - fdy 1 ^ = 0, 
(f)dy°W 3 - jdy x W 2 = 0. 



Remember that W3 is a scalar and W 2 is a pseudo-scalar, the above equations strongly 
suggest that these quantities are constant. If one further inserts the known result W 2 = 
into above equations, then § dy°W 3 = § dy x W 3 = 0. 
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The conclusion from the study of another set of MM equation and Ward identity is 
that Ws = is altogether a reasonable assumption. Indeed we believe that our argument 
presented above is close to a proof. The fact that both (B.3) and (B.5), which have 
been crucial for our consistency check, are valid only when x is not an intersection point, 
cautions us that W3 = may be not true at an intersection point. Luckily, for our solution 
of the Wilson loop in sect. 3, we do not need touch intersection points. 
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